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Abstract 

We construct analytical and numerical vortex solutions for an extended Skyrme- 
Faddeev model in a (3 + 1) dimensional Minkowski space-time. The extension is 
obtained by adding to the Lagrangian a quartic term, which is the square of the 
kinetic term, and a potential which breaks the 50(3) symmetry down to 50(2). 
The construction makes use of an ansatz, invariant under the joint action of the 
internal 50(2) and three commuting U{1) subgroups of the Poincare group, and 
which reduces the equations of motion to an ODE for a profile function depending 
on the distance to the x^-axis. The vortices have finite energy per unit length, and 
have waves propagating along them with the speed of light. The analytical vortices 
are obtained for special choice of potentials, and the numerical ones are constructed 
using the Successive Over Relaxation method for more general potentials. The 
spectrum of solutions is analyzed in detail, specially its dependence upon special 
combinations of coupling constants. 
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1 Introduction 



The so-called Skyrme-Faddeev model was introduced in the seventies [T] as a generali- 
zation to (3 + 1) dimensions of the 0(3) non-linear sigma model in (2 + 1) dimensions 
[2]. The Skyrme term, quartic in derivatives of the field, balances the quadratic kinetic 
term and according to Derrick's theorem, allows the existence of stable solutions with 
non-trivial Hopf topological charges. Due to the highly non-linear character of the model 
and the lack of symmetries, the first soliton solutions were only constructed in the late 
nineties using numerical methods [3l HJ |5l [6] . Since then the interest in the model has 
increased considerably and it has found applications in many areas of physics due mainly 
to the knotted character of the solutions [7J. The numerical efforts in the construction 
of the solutions have improved our understanding of the properties of the model [8] and 
even the scattering of knotted solitons has been investigated [9]. One of the aspects of 
the model that has attracted considerable attention has been its connection with gauge 
theories. Faddeev and Niemi have conjectured that it might describe the low energy limit 
of the pure SU{2) Yang-Mills theory [10]. They based their argument on a decomposition 
of the physical degrees of freedom of the SU{2) connection, proposed in the eighties 
by Cho ^T^, and involving a triplet of scalar fields n taking values on the sphere S"^ 
{n^ = 1). The conjecture, which is quite controversial [l2], states that the low energy 
effective action of the SU{2) Yang-Mills theory is the Skyrme-Faddeev action, and the 
knotted solitons would describe glueballs or even vacuum configurations. The fact that 
the Skyrme-Faddeev model has an 0(3) symmetry, and so possesses Goldstone boson 
excitations, is one of the many difficulties facing the conjecture, and some modifications of 
it were in fact proposed ^13]. Any check of such type of conjectures is of course very difficult 
to perform since it must involve non-perturbative calculations in the strong coupling 
regime of the Yang-Mills theory. However, Gies [H] has calculated the Wilsonian one loop 
effective action for the pure SU (2) Yang-Mills theory assuming Cho's decomposition, and 
found that the Skyrme-Faddeev action is indeed part of it, but additional quartic terms in 
the derivatives of the triplet n are unavoidable. In fact, the first numerical Hopf solitons 
were first constructed for the Skyrme-Faddeev model modified by a quartic term [3] which 
is the square of the kinetic term. However, the soliton solutions in [3j were constructed 
for a sector of the theory where the signs of the coupling constants disagree with those 
indicated by Gies' calculations. The addition of quartic terms has the drawback of making 
the Lagrangian dependent on terms which are quartic in time derivatives and so the energy 
is not positive definite. However, as a quantum field theory the Skyrme-Faddeev model is 
not renormalizable by power counting and has to be considered as a low energy effective 
theory. In addition, under the Wilsonian renormalization group fiow the square of the 
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kinetic term is as unavoidable as the Skyrme quartic term. Therefore, it is quite natural 
to investigate the properties of the Skyrme-Faddeev model with such modifications. 

In this paper we consider an extended Skyrme-Faddeev model defined by the La- 
grangian 

£ = d^n ■d''n-\ {d^n A d^nf + ^ {d^n ■ d^nf - V (rig) (1.1) 

where is a triplet of real scalar fields taking values on the sphere 5*^, M is a coupling 
constant with dimension of (length)"^, and f3 are dimensionless coupling constants, 
and the potential is a functional of the third component of the triplet n. Note that 
the potential breaks the 0(3) symmetry of the original Skyrme-Faddeev down to 0(2), 
the group of rotations on the plane ni n2, and so eliminating two of the three Goldstone 
boson degrees of freedom. In this paper the main role of potential is to stabilize the vortex 
solutions. 

The first exact vortex solutions for the theory (II. ip were constructed in [15] for the 
case where the potential vanishes, and by exploring the integrability properties of a sub- 
model of (11. ip . In order to describe those exact vortex solutions it is better to perform 
the stereographic projection of the target space onto the plane parameterized by the 
complex scalar field u and related to n by 

n= {u + u*, -i {u-u*),\u p -1) / (1+ I M P) (1.2) 

It was shown in [iBj that the field configurations of the form 

u = u{z,y) u* = u*{z*,y) for /3e^ = l 1^ = (1.3) 

are exact solutions of (11. ip . where z = + isix^ and y = — e2x'^ , with Sa = ±1, 
a = 1, 2, and x'^,/i = 0,l,2,3, are the Cartesian coordinates of the Minkowski space-time. 
Despite the fact that (II. 3p constitutes a very large class of solutions, no finite energy 
solutions were found within it. If the dependence of the u field upon the variable y is in 
the form of phases like e*'^^, then one finds solutions with finite energy per unit of length 
along the x^-axis. The simplest solution is of the form u = e*'^^, with n integer, and it 
corresponds to a vortex parallel to the x^-axis and with waves traveling along it with the 
speed of light. More general solutions of the class (II. 3p were constructed in [16], including 
multi-vortices separated from each other and all parallel to the x^-axis. The ideas of 
[T5] were generalized to an extended Skyrme-Faddeev defined on the target space CP^ , 
possessing — 1 complex scalar fields Ui, and the class of solutions constructed is like 
(II. 3p . where the fields Wj's are arbitrary functions of z and y pLTj. Note that the solutions 
(II. 3p are not solutions of the original Skyrme-Faddeev model, since that corresponds to 
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(3 = 0, and (11 .Sp requires the condition /3 = 1. If one takes the hmit /3 — )■ and ^ — J- 
with keeping the product /3 constant and equal to unity, one observes that (11.11) reduces 
to the CP^ model (if ^ = 0). Therefore, the configurations (11.31) are also solutions of the 
four dimensional CP^ model. The ideas of [17] were then used to construct multi vortex 
solutions for the the four dimensional CP^ model [HI [19] . 

The static energy density ("Hstatic = — -C) associated to (11.11) is positive definite if 
> 0, > 0, > and (3 < 0. That is the sector explored in [3] and where Hopf 
soliton solutions were first constructed (for V = 0). In addition, that is also the sector 
explored in [20] but with additional terms involving second derivatives of the n field, 
and where Hopf solitons were also constructed. The static energy density of (II. ip is also 
positive definite for > if 

M^>0; e^<0; /3<0; (3e^>l (1.4) 

That is the sector that agrees with the signature of the terms in the one loop effective 
action calculated in [T3] and that we will consider in this paper. Static Hopf solitons were 
constructed in [211 [22] for the sector (ll.4p (with = 0) and their quantum excitations, 
including comparison with glueball spectrum, were considered in [23]. An interesting fea- 
ture of the Hopf solitons constructed in [21] is that they shrink in size and then disappear 
as /3e^ — )■ 1, which is exactly the point where the vortex solution of the class (II. 3p exists. 

The aim of the present paper is to investigate if vortex solutions for the model (II. ID 
continue to exist when the condition /3 = 1 is relaxed, and so if they co-exist with the 
Hopf solitons in [21] . We also aim at the study of their properties and stability. The idea 
is to keep the solutions as close to those of the class (I1.3P as possible. In order to do that, 
we follow the ideas of [21] and implement an ansatz based on the 0(2) internal symmetry 
given by the transformations u e*" u, together with three commuting transformations of 
the Poincare group given by rotations on the plane x^, and translations in the directions 
x^ and x^ . We impose the field configurations to be invariant under the diagonal subgroups 
of the tensor product of the internal 0(2) group with each one of the three commuting 
one parameter subgroups of the Poincare group. The resulting ansatz is given by 

« = /(p) e'("^+^ (1.5) 

where n is an integer for u to be single valued, and A, k are real dimensionless parameters, 
and where we have used the dimensionless polar coordinates (p, ip, z, r), defined by 

x^ = ct = TqT x^ = rQpcosip x^ = ropsin(^ x^ = z (1.6) 

and where we have introduced a length scale tq given by 

4 , . 



which is positive since we are deahng with < (see (11. 4p ). Note that when A = ±k 
and / ~ p^", the configurations (ll.Sp are of the type (ll.3p . The ansatz (II .Sp is in fact 
a generahzation to (3 + 1) dimensions of the ansatz used in the Baby Skyrme models 

ESI ESI. 



The types of potential we will consider in this part are of the form 

V (ns) = Y (1 + nsf-^ (1 - ^3)'+^ (1.8) 

where is a non- vanishing integer, fi a real coupling constant. It is interesting to note 
that when the integer n of (II. 5p has the same modulus as N of (11.80 . one obtains analytical 
solutions of the form 



with e = ±1, and where a =\ N 



(-e2)(/3e^-l)M4 



ll/4 

. Such exact solutions are valid for 



all values of the coupling constants. In particular, for the case /3 = 0, (II. 9p are exact 
solutions of the theory 

C = M^ d^n -d^n-^ {d^n A d,nf - ^ (1 + ng)'"* (1 - ^3)'+^ (1.10) 
which is the proper Skyrme-Faddeev model in the presence of a potential. In this case we 



have a =\ N 



1 1/4 

and /i^ > 0. 



In some cases we have been unable to find a numerical solution which is expected in the 
analytical set, i.e. (II. 9p . Reasons for this are presented below. Apart from those cases, we 
have checked numerically the existence of the above solution. The numerical simulations 
are performed using a standard technique for a differential equation, the Successive Over 
Relaxation (SOR) method. In order to further confirm the accuracy and correctness of the 
SOR code, some of the results were reproduced by an independent code using Newton's 
method, giving typical differences of the order of less than 10~^. 

The paper is organized as follows. In the next section we briefly describe the extended 
Skyrme-Faddeev model. The equations of motion are also introduced in Sec. El We 
discuss the Hamiltonian density of the model in Sec. [3l The method and the solutions of 
the integrable sector of the present model are discussed in Sec. HI In Sec. [5l we show the 
numerical solutions. A brief summary is presented in Sec. [61 



2 The model 

In terms of the complex scalar field u introduced in (II. 2p the Lagrangian (II. ip becomes 
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£ = 4 . + — 

(1+ I u e2 



[1+ I U 



(1+ I n 



-V{\u\') (2.1) 



where we have used the fact that is a functional of | u p only, and so is the potential. 
The Euler-Lagrange equations following from (12.11) . or (II. ip . reads 



(1+ I u p) 5^/C^_2u*/C^a'^M = (1+ I u \y V 



(2.2) 



where V 



dV 



a\u\ 



^1 



and 



r=A#2^ , 4 [(.9.^g'-^).9^»* + (/3e^ - 1) (^.n^'^n'-) g^n] 
"""^ (1+UP)2 



(2.3) 



We point out that the theory (12.11) possesses an integrable sector defined by the condition 



{d^^uf = 



(2.4) 



Such condition was first discovered in the context of the CP^ model using the generalized 
zero curvature condition for integrable theories in any dimension [27], and then applied 
to many models with target space being the sphere S*^, or CP^ (see [28] for a review). It 
leads to an infinite number of local conserved currents. Indeed, (12.41) together with the 
equations of motion (12. 2p imply the conservation of the infinity of currents given by 



J. 



G 



/C, 



5G_ 

5u 



r* — 
^ 5u* 



(2.5) 



where G is any functional of | m p only. For the case where the potential vanishes, the 
set of conserved currents is considerably enlarged since G can be an arbitrary functional 
of u and u* , but not of their derivatives. If in addition to the condition (12. 4p one takes 
y = and /3e^ = 1, then the equations of motion reduce to d'^u = 0. It is in that 
integrable sector that the solutions (II. 3p lie, and were studied in [15j. For theories defined 
by Lagrangians which are functionals of the Skyrme term only (pullback of the area form 
of the sphere) the currents of the form (12. 5p are Noether currents associated to the area 
preserving diffeomorphisms of 5*^ [29] . It is possible to define conditions weaker than (12. 4p 
that lead to integrable theories associated to Abelian subgroups of the group of the area 
preserving diffeomorphisms [30] . 

Substituting the ansatz (II. 5p into the equations of motion (12.20 we get an ordinary 
differential equation for the profile function / as 
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(1 - P) 



s 



A- 



2\2 



4M2 



(1 + /^) 



dV 

d \ u\ 



(2.6) 
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where the primes denote derivatives w.r.t. p, and where 

t2 



n 



S = l+I3e 



R=l + 



2 f 



With the choice of potential given in fll.Sp we get that 



4M2 



d I u 



M2 



(2-^) (/^)(^--)-(2 + |) 



We look for solutions satisfying the following boundary conditions 



n — )• 



(0,0,-1) forp^O 



(0,0, +1) forp^cx) 
which imply that the profile function should satisfy 



/ for p 



and 



/ — !■ oo for p — )■ oo 



Let us then assume the following behavior of the profile function 



/ 



ao p'^^ (l + «! p + 0:2 p^ . . . ) for p — )■ 
/3o/^ (l + /3ii + /32^...) forp-^00 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



where Sj > 0, i = 1, 2. Substituting that into the equation (12.61) one gets the behavior for 
small p implies that 



si 



n 



(2.12) 



where n is the integer in the ansatz fll.51) . The behavior of fl2.6p for large p implies that 
the relation between v? and s| depends upon the form of the potential, and 
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M2 



for = -1 

for = -2 
for all other A^ 



(2.13) 
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where N is the integer appearing in the potential (11. 8p . and A, k are the parameters of 
the ansatz (II. 5p . Therefore, except for the cases = — 1 and N = —2, the waves along 
the vortex have to travel with the speed of light since the dependency upon and x° 
has to be of the form ± ct. For the dimensionfull constants L := TqA, and K := r^k, 
the velocity is defined as 



f , ^\ <c for = -1 
L I > c for AT = -2 



(2.14) 



Therefore, the mode is tachyonic for A^ = —2. A^ = — 1 is not tachyonic, but the energy 
diverges from the boundary behavior of the potential. In the following analysis, we will 
concentrate on the analysis for A^ ^ 1 (thus = /c^). 



3 The Energy 



The Hamiltonian density associated to (12. ip is not positive definite due to the quartic 
terms in time derivatives. We shall arrange the Legendre transform of each term in (12. ip 
to make explicit such non positive contributions, and write the Hamiltonian density as 
(see [32] for details) 



n = 4:M^ 



24 



ii P +Vu ■ Vu* 



(1- 



24 

^2 



Vu{ (Vu*^^ 



(1- 



u 



We' -I) 


1 M P +1 Vm ■ Vu* 




Vu ■ Vu*— u p 


e2 


(1+1 





(3.1) 



V{\u\') 



where ii denotes the x'^-derivative of u, and Vu its spatial gradient, and where we have 
denoted 



Vu 



1 + 



(3.2) 



with F > and < $ < 27r, being functions of the space-time coordinates. Note that H 
given in (13. ip is positive definite for static configurations and for the range of parameters 
given in (II. 4p . 

Using the ansatz (II. 5p and the coordinates (II. 6p we get u = iku/rQ. The metric on 
the spatial sub-manifold is given by ds' = {dp' + p' dip' + dz'), and so 



Vu 



u 



(^]_ - A^) 



Vu ■ Vu* 



{n^ + A^) 



(3.3) 
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where 



± 



f 



In addition one gets that , " . ^ = I + F e* * 

, 2 



fc2 

■(n_A2) 



(3.4) 

and since it is real it follows that 



$ = or TT. Therefore, (§)" - = (^]_ - - 3^) (^]_ + P - A^) /3 (^]_ - A^)^ So, 
the Hamiltonian density (13. ip becomes 



r, 







+ A^ + A;2) 



+ 2 , 1-1 (r]_ - A^ - 3fc2) (r]_ + - A^) 



(3.5) 



4 The integrable sector 



It is interesting to note that (I2.6l) -( l2l9l) with a special choice of the potential (ll.Sp have 
an analytical solution for each topological charge n. In fact, solutions of the integrable 
equation (12. 4p also become the solutions of the present model. For X'^ = k"^, the solutions 
can be written of the form 



P\ ""pile nip+k{z+T) 

a) 



(4.1) 



where e = ±1, and a is a dimensionless constant to be fixed by the equations of motion. 
Substituting this into the equation (I2.6p - (l2.9p . we get 



We' - 1) 



2^ 

M2 



4n^ 



1 



a) 

2N^ -3 



"[(n-l)n 



2n~4 

a J 

2N-A 



(n + 1) 



<^4)(^ 



An—A- 



AN-A-i 



(4.2) 



The constant a determines the scale of the vortex and the equation is satisfied if ri = 
and 



a =\ n 



'M\(3e'-iy 


1/4 


= n 





4/i2 



1/4 



(4.3) 



Thus, for all possible values of /3e^ we have analytical solutions. All those solutions satisfy 
the condition (12.40 . Clearly the class of solutions contain the special solution at /3e^ = 1 
found previously in [15] if we take a proper limit of vanishing potential, i.e. /3e^ — )■ 1 
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and fx^ — )■ 0, with = constant. Also, apparently we have no solution at (3e^ ^ 1 

without any potential because the scale (14. 3p goes to infinity. Note that the case /3 = is 
particularly interesting since it corresponds to the proper Skyrme-Faddeev model (without 
the extra quartic term) in the presence of a potential. Therefore, the configurations (14. ip 



4At2 



1/4 

, and 



are exact solutions of the theory fll.lOp (for n = N), with a =\ N 
/i2 e2 > 0. 

As we mentioned in SeclH for the sector satisfying (12. 4p . the model possesses the 
infinite set of conserved currents (12. 5p . In particular, choosing a form of G = — 4i(l + 
one gets of the Noether current for the symmetry of an arbitrary angle a, i.e., 
u — )■ e'^M 

(1 + \u\^)'^ (1 + r^l ) 

For the solution (14.11) . we can evaluate the charge per unit length for the solution 

/r n 1 1 

dxidxiJo = -87TM^ka\o lin) + — -(/^e^ - 1) (4.5) 
L 6 J 

where J(n) = ^r(l + ^)r(l — ^), with F being the Euler's Gamma function. Here we 
used an integral formula [3l] 

x^'^'dx 1 ^j9^i^F(^)F(l + m-^) ^^^^ 



^.iaiic = 27r + — -(/3e^-l) (4.7) 



(p + gx^)'"+i up^'+^\qJ F(l + m) 

For the Hamiltonian (13. 5p . we perform the similar computations. As a result, we get 
the energy per unit length by the integration on the X1X2 plane. For n = 1, the energy of 
the static vortex is (in units of 4M^) 

47r 1 

and for n>2 they are 

Estatic = 27in + - - l)^(^) (4-8) 

3 

Note that the first term is proportional to the topological charge. The energy per unit of 

length of the time-dependent vortex diverges for n = 1, and for n ^ 2 we obtain 

Ovr 1 r 27r 1 

E^ave = iT^n + —-A^e^ - \)[n^ - \)l[n) + P n^aH[n) + — (/Je^ - \)n (4.9) 
6 L 3 J 



For the limit of /3e^ — )■ — )■ with keeping = ^/M^fie^^^^Y)Jrl^ finite, we 
obtain the energy per unit of length found previously in [15] . The energy monotonically 
grows as k"^ increases. Interestingly, the static vortex has a minimum at /3e^ = 1.0 and/or 
/i^ = 0.0 but for the time dependent vortex there is a minimum of the energy for fixed /3e^ 
and k"^ and finite /x^. The solutions are confirmed numerically in the subsequent section. 
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5 The numerical analysis 



Although the ansatz (11.51) is given in terms of the polar coordinates, for the numerical 
analysis it is more convenient to use a new radial coordinate j/, defined by p = sj^-j^- 
Accordingly, we adopt a function g called the profile function, instead of using /, i.e., 

/(P) = 



9{y) 



The equation (12. 6p can be promptly rewritten as 

2ya(l-j,)\ \g(l-g)" 

^^''i-jj'-V-ij^- 2(1 + 1 



d 

dy 



9{^-9) . 



?/2 M2 

where the primes at this time indicate derivatives w.r.t.y and where 



y 



S = l + I3^g{l - g) 



R = l + g{l-g) 



1-y 



i-y 



9{l-. 
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(5.1) 

(5.2) 
(5.3) 
(5.4) 



The energy in the unit of AM'^ per unit length for the time-dependent vortex can be 
estimated in terms of following four parts of integrals of the dimensionless Hamiltonian 
H := n/AM^ 



E = 2tt pdpH{p) = E2 + ^ + iPe^ - 1)eP + 
Jo 



(2) , ^ 
M2 



En 



(5.5) 
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in which the components are defined as 



2(1 y v^?(i 



X 



(^^-A^)^+n^-(^.') K^a-^)r (5.7) 



""io 2(l-,)4^ + % + "^Ul-^^)" 



X 



{(.-A^)1^..^.(^.')}(.(1-.))^ (5. 



£;o = 27r l'%g^^^{\-gf-^ (5.9) 

For the integrable sector, we should choose N = n. 

Generally speaking, vortex is an object in three spatial dimensions, thus we have 
explored solutions in three spatial dimensions of (11.11) . In the three spatial dimensions, 
the 4th order terms in the Lagrangian (including the Skyrme term) successfully avoid 
the non-existence theorem of static and finite energy solutions by Derrick. However, the 
equation (12.61) of the ansatz (11. 5p is the same as an equation of corresponding static two 
spatial dimensions. This means z component has no essential contribution to the stability. 
In fact, the Derrick's theorem for two spatial dimensions implies that the contribution to 
the energy per unit length from quartic terms and the potential must be equal, namely 

(/3e^ - l)i?f = ^i?o (5.10) 

Since the solutions of the integrable sector satisfy E^^'^ = 0, putting together with /3e^ = 1 
and fi"^ = 0, one can confirm that the solution without the potential found in [15] satisfies 
the above condition automatically. This fact indicates that the Derrick's argument to the 
energy per unit length also works well outside of the integrable sector. 

The definition of the scale parameter a given in (14.31) indicates the existence of the 
analytical solution for the same sign of /3e^ — 1 and /i^. In our previous study of Hopfions 
on the extended Skyrme-Faddeev model, we confirmed numerically the solutions exist 
only for /3e^ > 1 [21], so we begin our analysis with the case of /3e^ > 1. We shall give 
comments for the possibility of finding solution of /3e^ < 1 in the next subsection. 
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5.1 The solutions of the integrable sector 



The analytic profiles f l4.ip can be written in the coordinate y as 



tor n 



jfe)=r*+aV* (5-11) 

where a is determined via (14. 3p . Apparently (15. lip are solutions of (15. ip . Next, we will 
see whether the solutions appear or not when we numerically solve (15. ip without any 
constraint. Also, for the obtained solutions we will check the zero curvature condition 

Since (15. ip is an ordinary second order differential equation, of course there are several 
methods to investigate. However, it is easily noticed that the equation (15. ip may exhibit 
singular-like behavior at the boundary because of the term g{l — g) of the denominator. 
Once the computation contains a small numerical error, the equation quickly diverges. 
The numerical method which can safely solve such a difficulty is well-known, the SOR 
method. Essentially we have solved the following diffusion equation for a field g{y,t) 



dg ^9 9^9 

5^ = "^" 



9, 



(5.12) 



in which we employ (15. ip as A. Here u is called as a relaxation factor which is usually 
chosen a; = 1.0 ~ 2.0. After a huge number of iteration steps, the field is relaxed to the 
static one, i.e, g{y,t) — > g{y), which we are finding. 

The first case is n = 1. From (II. 8p . the explicit form of the potential is 

^=1 = ^(1-^3)' (5.13) 

In Fig{T]we present the numerical solution g{y) and the corresponding Hamiltonian density 
H[p) for n = 1. Figj2]is the energy per unit length and its components for several values of 
(ie^ with fixed /x. The function g{y) in FigH] perfectly agrees with the analytical solution 
(15. lip . We shall give a few comments for the components of the energy. For the integrable 
solution, the topological contribution of the energy, i.e., E2 should be a constant. Also, 
E^^ is exactly zero for the integrable solution. The value of the component E^^'^ in the 
numerical solution is not exactly zero, but compatible with zero within the numerical 
precision. Note that the plot seems to blow up for the vicinity of /3e^ = 1.0, but the value 
is still up to order ~ 10~^, so it is still negligible. This clearly means that our numerical 
solutions satisfy the zero curvature condition and thus belong to the integrable sector. 
These numerical errors are probably originated in the finite number of the mesh points. 
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Figure 1: The n = \ profile g[y) and the corresponding Hamiltonian density of the real 
space -ff(p) of = 0.0 for the constant rg/i^/M^ = 1.0. 

In a usual case, we used the number A^mesh = 1000. When we employ a larger number, 
the value of E2 should be converge to the constant, i.e., 2ii. 

For the n = 2, form of the potential is 

Vn=2 = ^{l + ns){l-n,f (5.14) 

thus the potential is zero at both the origin and the infinity. Figj3]is the profile function 
and the Hamiltonian density for n = 2. Again the numerical profile and the analytical 
one (15.111) coalesce. Contrary to the case of = 1, the density has annular shape . FigJH 
is the energy per unit length and its components for several values of /3e^ and fixed /i^. 
Again we confirmed that the value of the component e[^'' is regarded as zero within the 
numerical uncertainty. 

For n = 3, form of the potential is 

K=3 = ^(l + n3)^/^(l-n3)«/^ (5.15) 

thus again the potential is zero at both the origin and the infinity. FigjS] is the profile 
function and the Hamiltonian density for n = 3. As is easily seen the radius of the annulus 
is larger than that of n = 2. FigJH]is the energy per unit length and its components for 
several values of /3e^ and fixed /x^. In this case, we face a numerical difficulty. During the 
computation by the SOR method, the solution g tends to oscillate around the true value 
and sometimes it accidentally goes below zero at the vicinity of the origin, and then the 
computation fails because of the term (yi^+^Z" in (15. ip . In order to avoid it, we employ a 
finer mesh, i.e., the number is at least A^mesh = 3000. 
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;ure 2: The static energy and its components corresponding to the solutions of FigJH 
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Figure 3: The n = 2 profile g{y) and the corresponding Hamiltonian density of the real 
space H{p) of fc^ = 0.0 for the constant rg/i^/M^ = 1.0. 

Until now, we have examined in the case of /3e^ > 1. The formalism leading to ( 14. 3p 
suggests that the choice /3e^ < 1 and < might also be possible and the scale is now 
defined as 

1/4 

for I3e^ < 1, /i^ < (5.16) 

The result is plotted in FigJTl However, existence of such solution seems dubious; the 
energy turns negative at a critical value of /3e^ thus the solution has no energy lower 
bound. Also, numerically the change of sign of the potential in the equation of motion 
quickly breaks the computation. 

5.2 The solutions outside of the integrable sector 

Although we have obtained the analytical solutions for a special form of the potential fll.Sp . 
we have many options for choice of the potential. The potentials which we employed in 
this paper essentially belong to a class of the generalized Baby-Skyrmion (BS) potential 
formally written as 

V^3s = y<, < = (l + n3r(l-n3r (5.17) 

For a = 0, the potential is called (the class of) the o/d-BS potential which was introduced 
in [33], while a 7^ is (the class of) the netoBS potential [M]. Our case (11. 8p corresponds 
to a = 2 — 2/n, 7 = 2 + 2/n. Note that the possible choice of the potential is certainly 



\n\ 



15 



30 
28' 
26' 
24- 
22- 
,^ 20 
18' 
16 
14' 
12 



1 ' r~ 

1 2 



~i 1 r 

9 10 



J3e- 




5ure 4: The static energy and its components corresponding to the solutions of Figj3l 
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Figure 5: The n = 3 profile g[y) and the corresponding Hamihonian density of the real 
space ii{p) of fc^ = 0.0. for the constant jM"^ = 1.0. 
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Table 1: The analysis of the limiting behavior of the solutions at both y = 0,1 for several 
choice of the potential, and O (^) which indicates that there exist (no) solutions. 

restricted by the analysis of the limiting behavior (12. lip for several potentials. The results 
are summarized in Table [H 

We can obtain many numerical solutions for the several types of the potentials. We 
show the result of n = 2 for the potentials f 2, "^4, "^^8/3; of course these are not of the form 
of the analytical solution. FigJH] presents the energies and the component E^^ for these 
potentials. For n = 2, the old-BS potentials give higher total energy than the new one. 
This indicates that the same class of potentials gives the similar energy and then, for 
n = 2 the energy of the new type potential v^j^ is closest to the integrable sector, which 
is also plotted in Fig. [8] for reference. 
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;ure 6: The static energy and its components corresponding to the solutions of FigJS] 
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Figure 7: The n = 1 profile functions and tlie energy density for the case of /3e^ < 1, for 
the constant Tq^u^/M^ = —1.0. 

6 Summary 

We have studied vortex solutions of the extended Skyrme-Faddeev model especially for 
the outside of the integrable constraint f3e^ = 1. In order to find the solutions, we 
introduced potentials of the extension of the Baby-Skyrmion type. We found several 
analytical solutions of the model. We also confirmed the existence of the solutions in 
terms of the numerical analysis. By using the standard SOR method, we obtained the 
axially symmetric solutions for charge up to n = 3 with several form of the potential for 
various value of the model parameters. 

In this work, we imposed the axial symmetry to the solution ansatz. However, solu- 
tions with lower symmetry, such as Z2-symmetry were found by a numerical simulation 
for the Baby-Skyrme model [35]. It would be interesting to investigate whether such de- 
formed solutions appear in the extended Skyrme-Faddeev model. Furthermore, full 3D 
simulations of the model will certainly clarify the detailed structure of the vortices. The 
analysis implementing these issues will be discussed in a forthcoming paper. 
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